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ABSTRACT 


This  thesis  discusses  a  procedure  that  optimizes  the  signal-to- 
noise  ratio  (SNR)  detected  by  a  linear  array  in  a  directional  ambient 
noise  field.  The  SNR  can  be  optimized  by  minimizing  the  ambient 
noise  detected  by  the  array.  For  a  given  target  location,  each 
possible  heading  of  the  array  centers  the  ambiguous  beam  of  the 
array  at  a  different  true  bearing.  Therefore,  each  heading  of  the 
array  will  receive  a  different  ambient  noise  level.  An  optimal 
heading  can  be  obtained  which  maximizes  the  SNR  received  by  the 
array.  For  all  possible  headings  the  beam  pattern  of  the  array  must 
be  determined  and  combined  with  the  ambient  noise  field  to 
determine  the  noise  level  detected.  This  thesis  discusses  the  theory 
involved  in  calculating  the  ambient  noise  levels  detected  for  each 
heading  and  provides  a  computer  program  which  performs  the 
calculations  for  a  particular  array. 
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1.  INTRODUCTION 


The  collapse  of  the  Soviet  Union  and  the  distribution  of  diesel 
submarines  to  many  third  world  countries  have  significantly  changed 
the  Antisubmarine  Warfare  (ASW)  problem  for  the  United  States  and 
its  allies.  The  open  ocean,  deep  water  ASW  scenario  versus  Soviet 
nuclear  powered  submarines  has  been  replaced  by  what  most 
experts  agree  to  be  a  shallow  water  scenario  against  quiet  diesel 
powered  submarines.  Additionally,  improvements  in  submarine 
sound  silencing  throughout  the  world  have  drastically  reduced 
passive  acoustic  detection  ranges.  The  change  in  the  ASW  problem 
and  the  improvements  in  sound  silencing  on  submarines  has  forced 
the  United  States  to  investigate  new  methods  to  increase  passive 
detection  ranges  while  maintaining  the  tactical  advantage. 

The  success  or  failure  of  a  tactical  towed  array  to  detect  a 
target  of  interest  often  depends  on  the  signal- to-noise  ratio  (SNR) 
received  at  the  sensors  making  up  the  array.  Any  opportunity  to 
increase  the  SNR  should  be  taken  advantage  of  to  increase  the 
probability  of  detection.  One  viable  means  of  increasing  the  SNR  is  to 
reduce  the  ambient  noise  detected  by  the  array.  The  ambient  noise 
detected  by  the  array  can  often  be  reduced  by  exploiting  the 
ambient  noise  field  directionality  by  orienting  the  array  in  the 
optimal  direction.  For  a  suspected  target  location,  each  array  heading 
centers  the  ambiguous  beam  (See  Section  IVA)  at  a  different  true 
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bearing.  Therefore,  each  heading  will  receive  a  different  ambient 
noise  level.  By  combining  the  beam  pattern  of  the  array  with  the 
ambient  noise  field  for  all  possible  array  headings,  an  optimal 
heading  can  be  obtained  which  minimizes  the  noise  level  detected  by 
the  array. 

The  majority  of  environmental  models  existing  in  the  fleet 
today  assume  an  omnidirectional  ambient  noise  field.  In  many  areas 
of  the  ocean  the  omnidirectional  model  is  sufficient,  while  in  other 
areas,  including  those  close  to  shore  and  shipping  lanes,  the  ambient 
noise  level  varies  significantly  with  true  bearing  in  the  frequency 
spectrum  dominated  by  distant  shipping.  The  directionality  of 
ambient  noise  can  be  used  to  the  tactical  advantage  of  the  operator 
by  orienting  the  linear  array  in  a  manner  which  aligns  the 
ambiguous  beam  so  that  the  total  ambient  noise  received  on  beams 
directed  towards  the  suspected  location  of  the  target  is  minimized. 
Linear  array  beam  pattern  calculations  and  the  interaction  of  these 
patterns  with  the  ambient  noise  field  are  discussed  in  detail  in  this 
thesis.  Additionally,  a  computer  program  is  provided  which 
performs  all  the  necessary  calculations  for  a  particular  array.  The 
program  is  limited  to  a  single  frequency  and  the  target's  suspected 
location  must  be  known  prior  to  running  the  program.  The  program 
is  provided  to  demonstrating  the  feasibility  of  constructing  a 
program  for  each  tactical  towed  array  in  the  Navy's  arsenal. 
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II.  HORIZONTAL  DIRECTIONALITY  OF  AMBIENT  NOISE 


A.  SOURCES  OF  DIRECTIONAL  AMBIENT  NOISE 

The  horizontal  directionality  of  ambient  noise  is  caused  by 
numerous  sources  including  rain  storms,  biologies,  and  distant 
shipping.  Most  of  these  sources  have  short  propagation  ranges  and 
rarely  affect  the  directionality  of  ambient  noise  detected  by  passive 
acoustic  sensors  except  for  distant  shipping.  The  frequency  range  in 
which  distant  shipping  is  the  principle  source,  and  therefore  subject 
to  high  levels  of  directionality,  is  from  about  50  to  500  Hz.  [Ref.  1]. 
Within  this  frequency  range  the  ambient  noise  can  be  dominated  by 
noise  which  can  travel  for  distances  of  1000  miles  or  more  from  the 
location  of  the  sensor  by  way  of  refracted  raypaths  in  the  deep 
sound  channel.  In  this  case,  the  ambient  noise  often  has  a  distinct 
directional  character  with  noise  levels  varying  between  true 
bearings.  It  is  also  within  this  frequency  range  that  many  passive 
acoustic  sensors  can  detect  enemy  submarine  signatures  including 
blade  rate  information. 

There  are  two  major  mechanisms  which  account  for  the  sound 
of  distant  shipping  being  trapped  in  the  deep  sound  channel  and 
therefore  contributing  to  the  directional  character  of  the  ambient 
noise  field.  The  first  mechanism  is  the  tendency  of  the  deep  sound 
channel  to  approach  the  surface  as  the  latitude  increases  [Ref.  2].  As 
can  be  seen  in  Figure  1 ,  the  depth  of  the  deep  sound  channel  in  the 
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Pacific  Ocean  ascends  toward  the  surface  as  the  latitude  approaches 
the  poles.  This  phenomenon  occurs  in  the  Northern  Atlantic  and 
Pacific  where  shipping  traffic  tends  to  concentrate  when  following 
great  circle  paths  across  the  oceans.  The  large  amount  of  shipping 
combined  with  the  shallow  deep  sound  channel  in  these  areas  cause 
a  directional  noise  field  in  the  vicinity  of  a  passive  acoustic  array 
which  is  located  in  the  sound  channel.  This  phenomenon  accounts  for 
the  higher  levels  of  directional  ambient  noise  fields  in  the  higher 
latitudes. 


Figure  1.  Depth  of  SOFAR  axis  in  the  Pacific  Ocean.  Contour  interval 
200  meters  [Ref.  3]. 

The  second  mechanism  accounting  for  shipping  noise  being  trapped 
in  the  deep  sound  channel  is  the  "downslope  conversion 
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phenomenon"  [Ref.  2].  This  describes  the  process  of  bottom 
reflected/surface  reflected  raypaths  changing  to  refracted  raypaths 
as  the  sound  crosses  over  the  continental  shelf  into  the  deep  ocean 
basin  and  into  the  deep  sound  channel.  Dense  shipping  normally 
found  close  to  shore  causes  this  phenomenon  to  cause  significant 
directional  noise  towards  the  shore  when  operating  off  the 
continental  shelf.  Figure  2  displays  both  of  these  mechanisms 
occurring  at  the  location  marked  on  the  chart.  The  lobe  to  the 
northwest  can  be  attributed  to  the  high  levels  of  shipping  traffic 
located  south  of  the  Aleutian  Islands  traveling  along  the  great  circle 
track  between  the  Far  East  and  North  America.  The  lobe  towards  the 
northeast  may  be  caused  by  the  down-slope  conversion  of  the  sound 
of  shipping  on  the  continental  shelf  off  the  west  coast  of  the  United 
States  [Ref.  4].  As  the  ASW  mission  moves  from  an  open  ocean,  deep 
water  scenario  to  a  scenario  closer  in  to  shore,  the  down-slope 
conversion  phenomenon  will  become  a  more  important  factor 
contributing  to  the  directional  nature  of  the  ambient  noise  field. 

B.  OBTAINING  THE  DIRECTIONAL  NOISE  FIELD 

The  directionality  of  the  noise  field  can  be  obtained  from  noise 
modeling  based  on  historical  data  or  from  real  time  in-situ 
measurements.  The  in-situ  measurements  have  the  highest  levels  of 
confidence  of  being  valid  for  a  given  scenario. 

1.  Noise  Modeling 

The  amount  of  noise  generated  by  distant  shipping  is 
directly  proportional  to  the  amount  of  shipping  traffic  present. 
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Figure  2.  Horizontal  Directionality  of  Ambient  Noise.  Site  location 
marked  by  X  on  top.  Directional  patterns,  as  obtained  with  end-fire 
beam  of  line  array  (dashed)  and  with  broadside  beams  (solid)  on 
bottom  [Ref.  2]. 


Therefore,  the  directionality  of  the  ambient  noise  field  can  be 
estimated  by  the  amount  and  location  of  the  shipping  in  the  vicinity 
of  the  sensor.  Historical  data  are  available  on  the  shipping  density 
around  the  world  and  can  be  used  to  model  the  expected  directional 
noise  field  for  a  given  location  at  a  given  time  of  year.  The 
environmental  factors  affecting  sound  transmission  in  the  oceans 
must  be  modeled  as  well. 

2.  Noise  Directionality  Measurements 

The  most  reliable  way  to  obtain  the  directionality  of  the 
ambient  noise  field  is  to  measure  it  with  on  board  sensors.  This  can 
be  done  by  performing  any  number  of  different  towship  maneuvers 
that  will  allow  the  measurement  of  beam  noise  data  on  three 
different  array  headings.  In  order  to  resolve  the  ambiguities,  each 
heading  must  be  at  least  forty  degrees  apart.  Wagstaff  s  Iterative 
Technique  [Ref.  5]  is  an  ambiguity  resolution  and  noise-field 
estimation  algorithm  which  can  be  used  to  perform  in-situ 
calculations  of  the  directional  noise  field. 
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III.  TOWED  ARRAY  BEAM  PATTERN 


A.  BEAMWIDTH 

All  equations  discussed  below  are  taken  from  Lawrence  J. 
Ziomek's  book  Underwater  Acoustics  [Ref.  6]. 

1.  Beamwidth  at  Arbitrary  Steering  Angle 

When  a  given  towed  array  beam  pattern  is  steered  from 
broadside  towards  end-fire,  the  beamwidth  increases.  A  minimum 
beamwidth  value  occurs  at  the  broadside  beam  providing  the  finest 
bearing  resolution.  To  compute  the  beamwidth,  first  assume  that  the 
array  is  linear  and  is  oriented  along  the  X  axis.  Next,  let  us  calculate 
the  beamwidth  for  the  horizontal  beam  pattern.  As  seen  in  Figure  3, 
'P’  is  the  steering  angle  in  degrees,  measured  counter  clockwise  from 
the  positive  X  axis.  Also,  let  ‘P+  and  represent  the  angles 
corresponding  to  the  3-dB  down  points  of  the  horizontal  beam 
pattern.  Note  that  Figure  3  shows  both  of  the  half-beamwidth  angles 
to  be  equal  to  A*P/2.  Although  this  is  not  true  in  general,  the 
resulting  equation  for  A'P  will  still  give  accurate  results. 

In  the  XY  plane  the  direction  cosine  u  with  respect  to  the 
X  axis  is  given  by 

u  =  cos('P),  (1) 

and  wiil  be  used  for  calculation  purposes.  Therefore,  let 

u+=cos(v+)>0  (2) 

u- =  cos('P-)  >  0  (3) 

where 
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and 


u+>  u-, 

=  »[''  -  (A^P/Z), 


(4) 

(5) 


iP-  =  »P'  +  (A»P/2). 


(6) 


Figure  3.  Representation  of  the  beam  steering  angie  'P'  and  the 
beamwidth  A'P  [Ref.  6]. 


Thus, 

since 


Au  =  u+  -  u-  =  2sin('P')sin(A'P/2)  >  0 


(7) 


cos(a-p)  -  cos(a+p)  =  2sin(a)sin(p)  (8) 


giving 


A'P=  2sin"l[Au/2sin('P')],  (9) 

where  a'P  is  the  beamwidth  in  degrees  of  the  horizontal  beam 
pattern  centered  at  the  angle  'P'.  Note  that  'P'  =  %/l  corresponds  to 
the  broadside  beam  pattern,  and  that  'P'  =  0  or  n:  corresponds  to  the 
end-fire  beam  pattern  where  Au/2sin('P')  =  <»  and  Eq.  (9)  is  not  valid. 
Equation  (9)  is  only  valid  when  lAu/2sin('P')l  ^  1.  Therefore,  a 
different  approach  must  be  used  to  determine  the  beamwidth  at 
end-fire  geometry. 
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2.  Beamwidth  at  End-Fire  Geometry 

Restricting  ourselves  to  the  horizontal  beam  pattern, 
consider  the  case  when  the  steering  angle,  ‘F’,  is  zero.  Letting 


u.  =  cos(‘P-)  >  0 

(10) 

and  substituting  *P'=0  into  Eq.  (6)  yields 

=  A'P/2. 

(11) 

For  this  case, 

u'  =  cos  ('P')  =  cos  (0)  =  1 

(12) 

and 

u-  =  1  -  (Au/2). 

(13) 

Combining  Eq.  (13)  with  Eqs.  (3)  and  (11), 

cos(A'P/2)  =  1  -  (Au/2) 

(14) 

or 

Avp  =  2cos-l(l  -  (Au/2)) 

(15) 

where  A'F  is  the  beamwidth  in  degrees  of  the  endfire  beams. 

B.  BEAMPATTERN 

Consider  a  linear  array  composed  of  N  (odd)  identical,  equally 
spaced,  amplitude  weighted  elements.  For  the  purposes  of  this 
thesis,  only  rectangular  amplitude  weights  will  be  discussed.  It  is 
important  to  note  that  the  theory  pertains  to  all  types  of  amplitude 
weighting.  However,  by  limiting  the  discussion  to  rectangular 
amplitude  weights,  a  closed  form  expression  can  be  obtained  for  the 
directivity  function  of  the  array. 
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1.  Directivity  Function 

The  directivity  function  of  a  iinear  array  can  be 
expressed  as  a  function  of  the  frequency  and  the  appropriate 
direction  cosine,  which  in  our  case  is  direction  cosine  u.  Figure  4 
shows  a  typical  normalized  beam  pattern  steered  in  the  direction  u  = 
u'  in  direction  cosine  space. 


ICts4(tu)l 


Figure  4.  Magnitude  of  a  typical  normalized  beam  pattern  steered  in 
the  direction  u  =  u'  in  direction  cosine  space  [Ref.  6]. 

The  closed  form  expression  for  the  normalized,  horizontal  directivity 
function  of  a  linear  array  with  rectangular  amplitude  weighting  is 

DN(f,fx)  =  sin(;tfxNd)/Nsin(3tfxd)  (16) 

where 

fx  =  uA,  (17) 
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N  is  the  number  of  elements  in  the  array,  d  is  the  inter-element 
spacing,  u  is  given  by  Eq.  (1),  and  xis  the  wavelength  of  the  sound 
waves.  Upon  steering  the  beam,  the  resulting  directivity  function, 
DN'(f,fx),  is  simply  a  translated  version  of  the  original  directivity 
function,  Diglf.fx),  that  is, 

DN*(f,fx)  =  DN(f,fx-fx’).  (18) 

Therefore,  when  steering  the  beam,  Eq.  (16)  can  still  be  used  by 
replacing  fx  with  fx-fx'  where: 

fx-fx'  =  (u-u')/X.  (19) 

2.  Beamwidth  in  Direction  Cosine  Space 

To  determine  the  3-dB  beamwidth,  Au,  in  direction  cosine 
space  for  the  normalized  directivity  function  given  by  Eq.  (16),  set 
sin(;iNdu4./X)/Nsin(ndu+A)  =  0.707  (20) 


or,  rewriting, 

where 

so  that 


sin(Nnx)  -  N(0.707)sin(nx)  =  0 


X  =  u+dA=  Aud/2X, 


Au  =  2Xx/d. 


(21) 

(22) 

(23) 


Upon  finding  the  solution  to  Eq.  (21)  for  x  and  substituting  the  result 
into  Eq,  (23),  the  3-dB  beamwidth  in  direction  cosine  space  is 
obtained. 
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IV.  OPTIMAL  LINEAR  ARRAY  HEADING 


The  towed  linear  array  is  one  of  the  most  valuable  sensors  on 
board  any  Naval  ASW  unit.  The  high  levels  of  directivity  obtained 
and  the  relatively  high  degree  of  mobility  combine  to  make  the 
towed  linear  array  an  effective  ASW  sensor.  These  two  factors  have 
also  led  to  the  common  practice  of  placing  the  narrow  broadside 
beam  in  the  suspected  direction  of  the  target.  This  tactic  has  worked 
well  in  the  fleet,  especially  when  the  ambient  noise  field  is 
omnidirectional.  However,  in  many  cases  where  the  ambient  noise 
field  is  directional,  this  tactic  can  lead  to  a  severe  degradation  of  the 
towed  array's  performance.  The  interaction  of  the  ambiguous  beam 
with  the  ambient  noise  field  can  lead  to  a  less  than  optimal  noise 
level  received  by  the  array. 

A.  AMBIGUOUS  BEAMS 

The  beamwidth  of  a  linear  array  beam  pattern  has  a  minimum 
value  at  broadside  and  increases  to  a  maximum  value  at  the  endfire 
beams.  This  conical  beam  pattern  inherently  presents  a  left-right 
ambiguity  problem  as  seen  in  Figure  5.  It  is  impossible  to  determine 
from  which  side  of  the  array  the  sound  is  coming  without 
maneuvering  the  array.  Each  beam  of  the  array  has  a  mirror  image 
on  the  opposite  side  of  the  array  commonly  called  the  ambiguous 
beam.  The  ambiguous  beam  receives  additional  noise  which,  in 
many  cases,  raises  the  total  noise  level  received  by  the  array  which 
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centers  the  main  lobe  of  the  beam  pattern  on  the  target.  It  is  often 
possible  to  minimize  the  ambient  noise  arriving  at  the  ambiguous 
beam  by  prudently  choosing  the  heading  that  the  array  is  oriented. 


nnt6en  of  Tow 


Figure  5.  Generic  Linear  Array  Beam  Pattern. 

B.  TOTAL  NOISE  LEVEL  CALCULATION 

When  numerous  individual  noise  sources  are  presented  to  a 
linear  array,  they  must  all  be  considered  in  determining  the  total 
noise  level.  Assuming  that  the  individual  noise  levels  are 
uncorrelated,  the  total  noise  level  may  be  calculated  using  the 
following  equation: 


TNL  =  10  loglZilO^^^i^l^b 


(24) 


where  NLi  are  the  individual  noise  levels 


noise  level  [Ref.  7]. 

C.  OPTIMAL  HEADING  DETERMINA 
Figure  6  displays  a  directional  amb 
^Pa  in  1  degree  increments.  The  target’s  si 
northwest. 


N 


Figure  6.  Illustration  of  the  interactioi 
ambient  noise  field. 


noise  averages  about  50  dB.  Using  Eq.  (24)  and  neglecting  detailed 
beam  pattern  calculations,  the  approximate  total  noise  level  received 
for  heading  1  is  about  62  dB.  Array  heading  2  is  to  the  southwest 
and  places  the  target  in  the  broadside  beam  (dashed  beam  pattern). 
The  ambiguous  beam  for  heading  2  is  located  to  the  southeast  where 
the  ambient  noise  level  averages  about  60  dB.  Again  using  Eq.  (24) 
and  neglecting  detailed  beam  pattern  calculations,  the  total  noise 
level  received  for  heading  2  is  about  67  dB.  Array  heading  1 
receives  a  total  noise  level  which  is  about  5  dB  less  than  that  of  array 
heading  2.  The  SNR  received  on  aiTay  heading  1  will  subsequently 
be  greater  than  the  SNR  received  on  the  broadside  array  heading  2 
and  the  probability  of  detection  will  likewise  be  proportionately 
higher.  Assuming  that  the  target's  signal  strength  does  not  vary 
significantly  over  the  time  it  would  take  to  maneuver  the  array, 
heading  1  should  be  chosen  over  the  broadside  heading  2. 

This  technique  can  be  used  to  determine  the  actual  noise  level 
detected  by  the  array,  including  detailed  calculations  for  the  beam 
pattern,  for  all  possible  headings.  The  total  noise  levels  detected  for 
each  heading  can  be  compared  with  the  noise  levels  for  all  the  other 
possible  headings  to  determine  the  heading  which  detects  the 
minimal  amount  of  noise.  Thus,  an  optimal  heading  can  be 
determined  for  the  orientation  of  the  array  which  maximizes  the  SNR 
detected  by  the  array  [Ref.  8]. 
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V.  PROGRAM  DESCRIPTION 


The  program  in  appendix  A  is  provided  to  demonstrate  the 
manner  in  which  the  optimal  heading  is  obtained  for  a  given 
directional  noise  field  and  for  a  suspected  target  location. 

A.  AMBIENT  NOISE  DATA 

The  program  prompts  the  user  to  input  the  ambient  noise  field 
by  one  of  two  means.  First,  the  user  can  input  the  name  of  a  file 
where  the  ambient  noise  data  is  stored.  This  would  correspond  to  a 
historical  type  of  ambient  noise  data  which  has  been  stored  in  the 
computer  at  a  previous  time.  The  second  means  of  inputting  the 
ambient  noise  data  is  for  the  user  to  input  the  data  directly  into  the 
computer.  This  process  is  somewhat  cumbersome  and  takes  about 
five  to  ten  minutes  depending  on  the  complexity  of  the  noise  field, 
but  the  data  is  copied  to  a  file  for  easy  access  at  a  later  time.  This 
corresponds  to  an  in-situ  ambient  noise  field  measurement. 

B.  LINEAR  ARRAY  CHARACTERISTICS 

The  program  simulates  a  linear  array  comprised  of  eleven 
omnidirectional  hydrophones  with  an  interelement  spacing  of  five 
meters.  The  interelement  spacing  allows  for  optimal  array 
performance  for  150  Hz  signals,  which  is  the  only  frequency  the 
program  is  capable  of  processing.  The  number  of  elements  generate 
a  broadside  beam  which  is  nine  degrees  wide  and  an  endfire  beam 
which  is  forty-six  degrees  wide.  The  characteristics  of  this  array  do 
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not  correspond  to  any  particular  array,  but  are  used  purely  to 
demonstrate  that  any  array  can  be  simulated  with  a  similar  program. 

C.  TARGET  INFORMATION 

The  only  information  that  is  required  concerning  the  target  is 
the  target's  suspected  true  bearing.  The  user  is  prompted  to  input 
this  information  which  the  program  uses  to  determine  the  relative 
bearings  to  the  target  for  all  possible  headings  that  the  array  may  be 
oriented.  Additionally,  the  program  determines  the  headings  which 
place  the  broadside  beam  in  the  direction  of  the  target.  The 
broadside  headings  are  used  as  a  reference  later  in  the  program, 

D.  BEAM  PATTERN  DETERMINATION 

For  each  possible  heading,  the  main  lobe  is  centered  at  the 
suspected  location  of  the  target.  The  program  also  determines  the 
bearing  at  which  the  ambiguous  beam  is  centered  for  each  heading. 
The  3-dB  beamwidths  of  the  beams  are  determined  using  either  Eq. 
(9)  or  Eq.  (15),  depending  on  the  orientation  of  the  target  to  the 
array.  The  program  also  determines  the  normalized  beam  pattern 
using  Eq.  (16).  The  true  bearings  which  the  main  lobe  and  the 
ambiguous  beam  cover  are  determined  for  use  in  the  noise 
calculation. 

E  NOISE  LEVEL  CALCULATION 

The  program  combines  the  beam  pattern  and  the  ambient  noise 
field  to  determine  the  ambient  noise  level  detected  for  each  possible 
heading.  The  ambient  noise  level  for  each  bearing  which  the  beam 
covers  is  converted  into  a  pressure,  multiplied  by  the  beam  pattern 
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factor,  then  summed  together  with  the  other  noise  levels  in  the  same 
beam  pattern  using  Eq.  (24).  After  all  the  pressures  are  summed  the 
total  noise  level  is  converted  back  into  a  dB  level  for  later  use.  For 
each  heading,  the  total  noise  level  detected  is  compared  to  the 
reference  broadside  heading  total  noise  level.  The  results  are  stored 
in  a  file  which  can  later  be  displayed  using  another  program. 

F.  PROGRAM  LIMITATIONS 

The  program  is  capable  of  simulating  one  linear  array  and  can 
only  process  noise  data  at  150  Hz.  Additionally,  only  rectangular 
windowing  is  used  when  determining  the  beam  pattern  of  the  array. 
Future  possible  work  could  provide  programs  which  can  simulate 
multiple  arrays  and  frequencies.  Other  types  of  windowing  could 
also  be  incorporated  into  the  program. 

G.  PROGRAM  EXAMPLE 

Figure  7  displays  an  ambient  noise  field  at  1 50  Hz  given  in  dB 
referenced  to  1  ^Pa.  The  noise  levels  are  given  in  one  degree 
increments  in  degrees  true.  The  higher  noise  levels  are  concentrated 
to  the  south  while  lower  noise  levels  are  found  to  the  north.  If  the 
suspected  location  of  the  target  is  to  the  north,  the  broadside 
headings  of  090  or  270  would  place  the  main  broadside  lobe  in  the 
direction  of  the  target.  While  these  courses  place  the  main  lobe  in 
the  direction  of  the  target,  the  ambiguous  beam  is  directed  towards 
the  south  in  the  direction  of  the  higher  ambient  noise  levels. 
Therefore,  courses  which  still  place  the  main  lobe  in  the  direction  of 
the  target,  but  place  the  ambiguous  beam  in  areas  of  lower  noise 
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AMBIENT  NOISE  PLOT 


Figure  7.  Directional  ambient  noise  field  at  150  Hz.  Levels  are  given 
in  dB  ref  1  ^Pa. 
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levels  may  have  lower  total  noise  levels  than  the  broadside  headings. 
Even  though  the  beamwidth  of  the  array  beam  pattern  is  greater 
than  the  broadside  beamwidth  and  more  individual  noise  levels  will 
be  detected,  the  total  noise  level  may  be  less  than  the  broadside 
heading  total  noise  level. 

The  optimal  heading  program  was  run  for  the  case  described 
above  and  the  results  are  displayed  in  Figure  8.  The  SNR  gain  for 
each  possible  heading  relative  to  the  broadside  heading  are 
displayed.  The  possible  headings  are  given  in  degrees  true  and  the 
SNR  gains  are  given  in  dB.  Note  that  for  the  broadside  headings  of 
090  and  270,  the  SNR  gain  is  zero  relative  to  broadside  which  is 
expected.  As  the  array  heading  is  changed  from  the  broadside 
heading  towards  the  north,  the  total  noise  level  received  by  the 
array  lowers  and  the  SNR  gain  relative  to  the  broadside  heading 
increases.  Note  that  when  the  array  is  oriented  about  30  to  60 
degrees  off  of  broadside  in  any  direction,  the  SNR  gain  is  from  5  to  6 
dB  relative  to  the  broadside  heading.  For  this  case  the  common 
practice  of  placing  the  broadside  heading  in  the  direction  of  the 
target  degrades  the  performance  of  the  array  by  about  6  dB. 
Additionally,  the  heading  to  orient  the  array  on  can  be  chosen  to 
comply  with  any  other  constraints  placed  on  the  ASW  unit. 

Additional  examples  using  various  ambient  noise  fields  are 
presented  in  Appendix  B  along  with  a  discussion  of  each  one. 
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SNR  GRIN  RELATIVE  TO  BROADSIDE  HERDING 


N 
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LEVELS  GIVEN  IN  OB  REF  IpPfl 
X  -  TARGET  TRUE  BEARING  360. 


Figure  8.  SNR  gain  for  all  true  headings  relative  to  the  broadside 
heading. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  optimal  array  heading  can,  in  many  cases,  significantly 
improve  the  SNR  for  a  linear  array  in  a  directional  noise  field.  Two 
major  applications  for  the  technique  occur  at  different  phases  of  a 
given  ASW  scenario.  The  first  application  occurs  when  planning  to 
search  an  area  for  a  target.  The  traditional  barrier  search  is  a  prime 
example  of  this.  The  target's  suspected  location  is  known  and  there 
is  plenty  of  time  to  consider  the  directional  ambient  noise  field  prior 
to  setting  up  the  barrier.  The  orientation  of  the  barrier  should  be 
based  upon  the  optimal  array  heading  calculation  for  the  given  noise 
field. 

The  second  major  application  occurs  after  a  target  has  been 
detected.  The  ASW  unit  should  place  itself  between  the  target  and 
the  direction  of  the  highest  ambient  noise  to  prevent 
counterdetection  from  the  target  submarine  while  optimizing  the 
unit's  detection  capabilities.  Additionally,  the  ASW  unit  should 
choose  the  courses  which  optimize  the  array's  SNR  based  on  the 
optimal  heading  calculation. 

B.  RECOMMENDATIONS 

Based  on  the  results  of  this  thesis,  the  following  actions  are 
recommended: 
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1.  A  computer  program,  such  as  in  Appendix  A,  should  be 
made  up  for  each  of  the  arrays  in  the  US  Navy's  arsenal  and 
distributed  to  the  applicable  ASW  units  capable  of  using  them. 

2.  A  fleetwide  ambient  noise  directionality  measurement 
procedure  should  be  made  up  and  distributed  to  the  fleet.  All  ASW 
unit  commanders  should  be  advised  to  the  practicality  of  measuring 
the  ambient  noise  field  and  using  the  results  to  optimize  the 
probability  of  detection. 

3.  A  directional  ambient  noise  data  base  should  be  established 
which  can  be  made  available  to  all  ASW  units. 
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APPENDIX  A 
FORTRAN  PROGRAMS 


This  appendix  provides  the  reader  with  the  program  "Optimal" 
which  is  compiled  using  the  FORTRAN  language.  The  program 
calculates  the  optimal  heading  to  search  for  a  target  in  a  directional 
ambient  noise  field.  The  program  is  described  in  detail  in  Chapter  V. 
Additionally,  two  plotting  programs  are  provided  for  display 
purposes.  The  first  program,  "Plot",  displays  the  directional  ambient 
noise  field  using  the  Disspla  routine  on  the  Naval  Postgraduate 
School's  mainframe  computer.  The  second  program,  "Plot2",  displays  ' 
the  results  of  the  Optimal  program  using  the  Disspla  routine  also. 
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FILEi  OPTIMAL  FORTRAN  A1 


DAVID  C.  MCDONNELL,  LT,  USN 

XXXXJfXXKXKXXXXXXXiiXXXXKXKXKXXXXKitXiOfXlIXXXidOOiKXKXXXIfXifilXKKXKXX 

XTHIS  PROGRAM  CALCULATES  THE  OPTIMAL  TOWED  ARRAY  HEADING  IN  X 
XA  DIRECTIONAL  NOISE  FIELD  ENVIRONMENT,  MINIMIZING  THE  X 

XAMBIENT  NOISE  DETECTED  BY  THE  SENSOR.  X 

XXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXKXXXXXKXXXXXXXXXXX 


VARIABLE  DECLARATION  SECTION 
CHARACTERXZO  FILEI , F1LE2, FILES 

REAL  X(3BO).Y(360).TOT,L(180),N1,N2,DN,HDO,DU,N,BH,OTHETA,DT1 
REAL  DT2, DT3, DTA, PSUM. PI ,ANS, LEVEL, BHEAD, PSY1P,PSY2P, F, FI 
REAL  G( 180 ) , PSYl , DPSYl , PSY2, DPSY2 
INTEGER  A,B,C,0.I,J.M.K,0,P,Q 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXX 

VARIABLE  DESCRIPTION  SECTION 

X  FlLEl,  FILE2,  AND  FILES  ARE  USER  FILES  TO  ACCESS  AND  STORE 
INFORMATION. 

XX(360)  IS  AN  ARRAY  WHICH  STORES  THE  AMBIENT  NOISE  DATA 
XY(360)  IS  A  DUMMY  ARRAY  HOLDING  THE  AMBIENT  NOISE  DATA 
XTGT  IS  THE  TARGETS  TRUE  BEARING 

KL(180)  IS  AN  ARRAY  WHICH  STORES  THE  AMBIENT  NOISE  LEVEL  FOR 
EACH  HEADING  BETWEEN  1  AND  180  TRUE 
XNI  IS  THE  NOISE  LEVEL  AT  THE  BEARING  OF  THE  TARGET 
XN2  IS  THE  NOISE  LEVEL  AT  THE  AMBIGUOUS  BEARING  OF  THE  TARGET 
XON  IS  THE  DIFFERENCE  IN  BEARING  BETWEEN  TARGET  AND  SENSOR 
XHOG  IS  THE  HEADING  OF  THE  SENSOR 

XDU  IS  A  CONSTANT  WHICH  AFFECTS  THE  BEAMWIDTH  OF  THE  ARRAY 
XN  IS  THE  NUMBER  OF  ELEMENTS  IN  THE  ARRAY 
XBW  IS  THE  BEAMWIDTH  OF  THE  ARRAY  IN  RADIANS 
KDTHETA  IS  THE  HALF-BEAMWIDTH  OF  THE  ARRAY  IN  DEGREES 
K0T1,0T2,DT3,DTA  ARE  THE  TRUE  BEARINGS  WITHIN  THE  MAIN  BEAN 
XPSUM  IS  THE  SUM  OF  THE  NOISE  PRESSURES 
XPIS3.1A1593 

XANS  3  ANSWER  TO  QUESTION 

XLEVEL  IS  THE  NOISE  LEVEL  IN  OB  OF  THE  AMBIENT  NOISE  REF  lUPA 
XBHEAD  IS  THE  BROADSIDE  HEADING  TO  THE  TARGET 
XPSYIP  AND  PSY2P  ARE  THE  ANGLES  IN  DEGREES  CCW  FROM  THE  SENSOR 
HEADING  TO  THE  TARGET  AND  THE  AMBIGUOUS  BEANS 
XF  IS  THE  BEAM  PATTERN  FACTOR 

XFI  IS  A  FACTOR  USED  IN  CALCULATING  THE  BEAM  PATTERN  FACTOR 
XGCIBO)  IS  THE  ARRAY  WHICH  STORES  THE  GAIN  FOR  EACH  HEADING 
XPSY1,PSY2, DPSYl, DPSY2  ARE  FACTORS  USED  TO  CALCULATE  THE  BEAM 
PATTERN  FACTORS 

XA,B,C.0,I, J,M  ARE  INTEGERS  USED  FOR  NOISE  FIELD  INPUT 
xK,0,P,q  ARE  INTEGERS  USED  FOR  THE  GAIN  CALCULATION 

PI=3. 141593 

XXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXKXXXXKXXXX 

OBTAIN  THE  AMBIENT  NOISE  DATA  FROM  THE  USER 

THE  USER  IS  GIVEN  A  CHOICE  TO  USE  EITHER  A  STORED  DATA  FIELD 
OR  TO  CREATE  HIS  OWN  DATA  FIELD 

WRITE(X,X)'PLEASE  CHOOSE  THE  NUMBER  CORRESPONDING  TO  THE  ' 
WRITE(X,X) 'MEANS  OF  ENTERING  THE  AMBIENT  NOISE  FIELD' 
WRITEIX.X)'  ' 

WRITE(X,X)'  1)  OBTAIN  DATA  FROM  AM  EXISTING  FILE' 

WRITE(X,X)'  2)  USER  INPUT  DATA' 

READ(X,x,ERR=5)  ANS 

IF  THE  USER  WANTS  TO  USE  A  STORED  DATA  FIELD 
IF(ANS.EQ.l)  THEN 

6  WRITE! K,x) 'ENTER  THE  NAME  OF  THE  FILE  WHERE  THE  NOISE' 

WRITE{x,x)'DATA  IS  STORED' 

READ<X, '(A10)',ERR=6)  FILEI 
OPEN(UNIT=l,  FILE=FILE1,  STATUS* 'OLD* > 


27 


OU  OM  UMN  a  M  K)  OUO  ^  «0  U^OUl^UUUUOUUUUOUOa  OOUOf^  uuoo 


FILE  I  OPTIHAL  FORTRAN  A1 
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DO  20  1^1,360 
READd.x)  X(I) 

CONTINUE 
CLOSE  (1) 

GO  TO  70 

IF  THE  USER  WANTS  TO  CREATE  HIS  OWN  DATA  FIELD 

ELSE  IF  (ANS.EQ.2)  THEN 

D=1 

M=1 

WRITE  («>«)  'ENTER  THE  FILE  WHERE  YOU  WISH  THE  NOISE' 
WRITE  (K.K)  'LEVELS  TO  BE  STORED' 

READ  (X. •(A10)',ERR=21)  FILE2 

OPEN  (UNIT=2,  FILE*FILE2,  STATUS* 'UNKNOWN' ) 


WRITE  (K,26)  D 

FORNATCENTER  THE  NOISE  LEVEL  BEARING' , IX, 13 . 3, IX, 'DEO  TRUE' 
WRITE  (X,X)>( ENTER  -99  TO  END  DATA  ENTRY  WHEN  REACH  BRO  360) 
READ  (K,X,ERR=25)  LEVEL 


) 


IF  (LEVEL .EQ. -99)  GO  TO  60 

WRITE  (K,X)>ENTER  THE  BEARINGS  WHICH  CORRESPOND  TO  THIS  LEVEL' 
WRITE  (X,30)  D 

FORMAT  ('STARTING  WITH' , IX, 13 . 3, IX, 'DEG  TRUE  (EX.  001,010)') 

READ  (x.X,ERR*29)  A,B 

C=B-A+1 

D*B+1 


THIS  DO  LOOP  WRITES  THE  NOISE  LEVELS  TO  THE  FILE 


DO  50  J*1,C 

WRITE  (2,90)  LEVEL 
FORMAT  (1X,F6.2) 
X(M)=LEVEL 
M»M+1 
CONTINUE 
GO  TO  25 


CL0SE(2) 

END  IF 

XXXXKXXXXXXXXXXXXXXKXXXXKXKXXXXKXXKXXXXXXXXKXXXKXXXXXKKXKKXXKX 

XXXKXXKKXKXXXXXKXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKKXX 


THIS  PROGRAM  USES  A  GENERIC  LINEAR  ARRAY  COMPOSED  OF  7 
ELEMENTS  SPACED  5  METERS  APART,  GIVING  THE  ARRAY  A  LENGTH  OF 
50  METERS.  THIS  SPACING  IS  SET  UP  FOR  A  FREQUENCY  OF  150  HZ. 
FOR  PURPOSES  OF  BEAMWIDTH  DETERMINATION  THE  FOLLOWING  EQUATION 
WAS  PREVIOUSLY  SOLVED  FOR  X. 

SIN(NXPIXX)-NX.7071XSIN(PIXX)=0 

RESULTING  WITH  X  *  .09091,  AND  DU=9XX  GIVING • 

DU*. 1616 
N*ll 


XXKXXXXXXXXXXXXXXXXXXKKXXMXXXXXXXXMXXKXXKXKXXXKKXXXXKKXKXXKXXX 

OBTAIN  THE  TARGET  INFORMATION  FROM  THE  USER 

WRITE(X,x)'EMTER  THE  SUSPECTED  TRUE  BRG  OF  THE  TARGET' 
READ(x,k,ERR*75)  TCT 
IF  (TGT.EQ.O)  TGT  *  360. 

DETERMINE  THE  BROADSIDE  HEADING  BETWEEN  1  AND.  180  FOR  THE 
GIVEN  TARGET  POSITION. 

BHEAD  *  TGT  -  90 
IF  (BHEAD. LE.O)  THEN 

BHEAD  =  TGT  ♦  90 
ELSE  IF  (BHEAD. GT. 180)  THEN 
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FILE:  OPTIHAL  FORTRAN  A1 


AHEAD  >  AHEAD  -  180 

XXKXMXitKXXKXKXXXXXXitKKXKXXXXKXXXXXXXXilXXXifKKKXXXXKKXXKliKXKKXXX 

OATAIN  NAME  FOR  FILE  WHERE  SNR  GAIN  DATA  IS  TO  AE  STORED 

WRITE  (X.X)'  ENTER  FILENAME  FOR  SNR  GAIN  DATA  STORAGE* 

READ  (X, '(AIO'.ERR^SO)  FILES 

OPEN  (UNIT  =  3,  FILE  =  FILE3.  STATUS^ *  UNKNOWN • ) 

WRITE  (3,X)  TGT 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXX 

FOR  EACH  POSSIALE  HEADING  THE  LOOP  CALCULATES  THE  GAIN 
RELATIVE  TO  AROADSIDE  THAT  THE  ARRAY  ACHIEVES 

HDG  =  1 

DO  200  K  =  1,180 
N1=TGT 

DO  90  0=1,360 
Y(Q)=X(q) 

CONTINUE 

0N=H0G  -  TGT 

THE  FOLLOWING  IF  STATEMENT  DETERMINES  THE  AMAIGUOUS  8EAN  FOR 
THE  CURRENT  HEADING 

IF  (DN.GT.-360.AND.DN.LE,-180)  THEN 
DH=DN+S60  •’ 

N2=HDG+DN 

ELSE  IF  (DN.GT. -180, AND.DN.lt. 180)  THEN 
N2=H0G+DN 

ELSE  IF  (DN.GE.180.AND.ON.LE.360)  THEN 

DN=S60-DN 

N2*HDG-DN 

END  IF 

IF  (N2.GT.360)  THEN 
N2=H2  -  360 

ELSE  IF  <N2.LT.O)  THEN 
N2=N2  +  360 
END  IF 

IF  (Nl.EQ.O)  Nl=360 
IF  (N2.Eq.O)  N2=360 

XXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXX 

DETERMINE  THE  SEARING  CCW  FROM  THE  ARRAY  FOR  THE  SEAM 
CONTAINING  THE  TARGET  AND  THE  AMAIGUOUS  SEAM.  THESE  RELATIVE 
SEARINGS  ARE  USED  IN  CALCULATING  THE  SEAMHIDTH  AND  THE  SEAM 
PATTERN  FACTOR  USED  SELOW 

IF  (HDG. GE. TGT)  THEN 
PSY1P=HDG-TGT 
PSY2P=360-PSY1P 
ELSE  IF  (HDG. LT. TGT)  THEN 
PSYlP=HDG-TGT+360. 

PSY2P=360-PSY1P 

END  IF 

XXXXXXXXXXXKXXXXXXXXXXXKXXKXXXXXXXXXXKKKXXXXXXXXXKXXXXXXXXXXXX 

DETERMINE  THE  SEAM  WIDTH  OF  THE  AEAM  CONTAINING  THE  TARGET  AS 
DELL  AS  THE  AMAIGUOUS  AEAM.  THE  3  DA  DOWN  POINT  IS  USED  TO 
DETERMINE  THE  SEAMHIDTH  OF  THE  AEAM  PATTERN. 
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C 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

92 

c 


c 

c 

94 

C 


FOR  ENDFIRE 

IF  (PSYIP.EQ.O. .OR.PSYIP.EQ.ISO.)  THEN 
BH  3  2xAC0S(l-(DU/2)) 

ELSE  IF  HOT  ENDFIRE 

ELSE  IF  (PSYIP.NE.O. .ANO.PSYIP.NE.ISO.)  THEN 

IF  (DU/(2»ABS(SIN(PSY1PXPI/’180))).LE.1)  THEN 

BM  >  2KASIN(OUy'(2XABS(SIN(PSYlP»PI/180)))) 
ELSE  IF  (DU/(2XABS(S1N(PSY1PSPI/180))).OT.1)  THEN 
BH  »  2XAC0S(l-(DU/2)) 

END  IF 

END  IF 

DETERMINE  THE  HALF  BEAHHIDTH  OF  THE  BEAM  IN  DEGREES. 
DTHETA  >  INT(BHX180/(PIX2)) 


XXXXXXXXXXXXXXXXKXXXXXXXXXXXXXKXXXXXXXXXXXXKXXXXXXKXXXKXXXXXXX 


CALCULATE  THE  TOTAL  NOISE  SEEN  BY  THE  ARRAY  ON  THIS  HEADING 
PSUN  «  10XX<X(INT(N1))/10)  *  10XX(X(1NT(N2))/10) 


THE  FOLLOHING  DO  LOOP  SUMS  THE  AMBIENT  NOISE  LEVELS  HITHIN  THE 
BEAMS  OF  THE  ARRAY.  THE  BEAM  PATTERN  FACTOR,  F,  IS  CALCULATED 
FOR  EACH  NOISE  LEVEL  HITHIN  THE  BEAM. 


DO  100  0«1, DTHETA 


NOISE  LEVELS  ARE  SET  TO  ZERO  AFTER  COUNTING  TO  PREVENT  THEM 
FROM  BEING  COUNTED  THICE. 

DTI  *  NI  +  0 

IF  (DTI. GT. 360)  DTI  »  DTI  -  360. 

PSYI»PSY1PXPI/180, 

DPSY1»<PSY1P+0)XP1/180. 

F1=(PI/2)X(C0S(PSY1)-C0S(DPSY1)) 

IF  (Fl.EQ.O)  THEN 
F=l. 

ELSE 

F>:SIN(NKF1)/(NXSIN(F1)) 

END  IF 

IF  (F.LT. 0.707)  GO  TO  92 

PSUM=PSUM+FxI0XX(Y(INT(0TI))/10) 

Y(INT(DT1))=0 


DT2  »  N2  ♦  0 

IF  (OT2.GT.360)  DT2  «  DT2  -  360. 

PSY2=PSY2PXPI/180. 
DPSY2=(PSY2P+0)XPI/180. 
Fl=<PI/'2)x(C0S(PSYZ)-C0S(DPSV2)) 
IF  (Fl.EQ.O)  THEN 
F»l. 

ELSE 

F«SIN(NKF1)/(NXSIN(F1)) 

END  IF 

1F(F.LT. 0.707)  GO  TO  94 

PSUM=PSUM>Fxioxx(Y(INT(DT2))/10) 

Y(INT(DT2))=0 


DT3  »  Nl  -  0 

IF  (DT3.LE.0)  DTS  >  0T3  *  360 

DPSY1*(PSY1P-0)xPI/I80. 
F1»(PI/2)X(C0S(PSY1 )-COS(DPSYl)) 
IF  (Fl.EQ.O)  THEN 
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F«I. 

ELSE 

F-SIN(NXF1)/(NXSIN(F1)) 

END  IF 

IF  <F.LT. 0.707)  GO  TO  96 
C 

PSUM=PSUM+F«10*K(Y(INT(DT3))/10) 

Y(IMTCDT3))=0 

C 

96  DT4  =  N2  -  0 

IF  (DT4.LE.0)  DT4  =  DT4  *  360 
C 

DPSY2»(PSY2P-0)«PI/180. 

FI  =  ( PI/2 ) *<C0S  <  PSY2 ) -cost  DPSY2 ) ) 

IF  (Fl.EQ.O)  THEN 
F=I. 

ELSE 

F=SIN(NKF1)/(NXSIN(F1)) 

END  IF 

IF  (F.LT. 0.707)  GO  TO  100 
C 

PSUM=PSUM+FX10XX( Y( INT(DT4 ) )/10) 

Y(INT(0T4))=0 

C 

100  CONTINUE 

C 

110  L(K)-10XLOG10(PSUM) 

c 

HDO  s  HDG  >1. 

200  CONTINUE 

C 

C  THE  FOLLOHING  DO  LOOP  CALCULATES  THE  GAIN  RELATIVE  TO 

C  BROADSIDE  FOR  EACH  HEADING. 

C 

DO  250  P  »  1,180 
C 

C  FOR  EACH  HEADING  THE  GAIN  RELATIVE  TO  THE  BROADSIDE  BEAM  IS 

C  CALCULATED  AND  PLACED  INTO  THE  ARRAY  G(P). 

C 

G(P)  »  L<INT(BHEAD))  -  L<P) 

C 

HRITE  (3,X)  G(P) 

250  CONTINUE 

C 

C  THE  FILE  HHERE  THE  GAIN  DATA  IS  STORED  IS  NON  CLOSED 

C 

CLOSE  (3) 

C 

END 
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FILE.  PLOT 


FORTRAN  A1 


DAVID  C.  MCDONNELL,  LT  USN  PLOOOOlO 
THIS  PROGRAM  USES  THE  DISSPLA  ROUTINE  TO  PLOT  THE  AMBIENT  NOISE  PL000020 
FIELD  BEING  USED.  THE  USER  MUST  IDENTIFY  THE  FILE  HHERE  THE  NOISEPL000030 
DATA  IS  STORED.  A  POLAR  PLOT  IS  USED  TO  DISPLAY  THE  NOISE  LEVEL  PLOOOOAO 
IN  DEGREES  TRUE.  PL000050 


IN  DEGREES  TRUE.  PL000050 

PL000060 

VARIABLE  DECLARATION  SECTION  PL000070 

PL000080 

INTEGER  I,J  PL000090 

REAL  X(360),Y(3<0}, 2(360), MIN, MAX, DR. DR1,DR2,0R3  PLOOOIOO 

CHARACTERK20  FILEO.FLNO  PLOOQllO 

PL000120 

OBTAIN  FILE  NAME  FROM  THE  USER  PLQ00130 

PLOOOIAO 

HRITE  (X,i()'£NTCR  FILENAME  AND  FILETYPE  HHERE  DATA  IS  STORED'  PL000150 

READ  (X.'(A20)')  FILE*  PLOOOIOO 

PL000170 

A  MUST  BE  APPENDED  TO  THE  FILENAME  FOR  THE  DISSPLA  ROUTINE  PL000180 

TO  RECOGNIZE  THE  FILE.  PLOOOIOO 

PL000200 

FLN*='/'//FILE$//'A'  PL000210 

PL000220 

OPEN  (UNIT»l.FILEsFLN$,STATUS*'OLD*)  PL0002S0 

PLOOOZOO 

THE  FOLLOHING  DO  LOOP  READS  THE  AMBIENT  NOISE  FROM  THE  FILE.  THE  PL0002S0 

DISSPLA  COMMAND  "POLAR”  PLOTS  THE  DATA  REFERENCED  TO  RADIANS  FROM  PL000260 


EAST.  THEREFORE,  THE  DEGREES  TRUE  ARE  TRANSFORMED  INTO  RADIANS 
FROM  EAST  FOR  PLOTTING  PURPOSES. 

INITIALIZE  VARIABLES 
MIN-IOO. 

MAX=0. 

DO  10  I»1.360 

READ  (I,x)  X(I) 

DETERMINE  THE  MINUHUH  AND  MAXIMUM  NOISE  LEVEL  FOR  PLOT  PURPOSES 

IF  (X(I).LT.MIN)  MIN=X(I) 

IF  (X(I).GT.MAX)  MAXsX(I) 

IF  (I.GT.0.AND.I.LT.90)  THEN 
J*90-I 

ELSE  IF  (I.GE.90)  THEN 
J*A50-I 
END  IF 

Y(J)=REAL<J) 

Z(J)=X<I) 

CONTINUE 

SET  MIM/MAX  EQUAL  TO  ♦/-lO  THAN  THE  ACTUAL  MIN/MAX  NOISE  LEVEL 
FOR  PLOTTING  PURPOSES. 

MIM=INT(MIN)-20. 

MAX^INTCMAXl-MO. 

DR=INT((MAX-MIN)/6.) 

DR1=MIN-»DR 
0R2  =  DR1-»DR 
ORSsDRZ^^OR 
DR4=DR3^DR 
DR5:DR4>DR 
MAXsDRS^DR 

DR=2KDR 

THE  PLOTTING  ROUTINE  USES  DISSPLA  COMMANDS  TO  PLOT  THE  DATA 
IN  DECREES  TRUE  HITH  THE  AMBIENT  NOISE  LEVELS  GIVEN  IN  DB  REF 
1  UPA. 

CLOSE(UNIT-l) 

CALL  SHERPA('MCT3NOIS','A',3) 


PL000270 

PL000280 

PL000290 

PL000300 

PL000310 

PL000320 

PLO003S0 

PL000540 

PL000350 

PL000360 

PL000370 

PLOOOS80 

PL000390 

PL000400 

PL000410 

PL000420 

PL000430 

PL000440 

PL000450 

PL000460 

PL000470 

PL000480 

PL000490 

PL000500 

PL000510 

PLOOOS20 

PL000530 

PL000S40 

PL000550 

PLO00S60 

PL000570 

PL000580 

PL000590 

PL000600 

PL0006IO 

PL000620 

PL000630 

PL000640 

PL000650 

PL000660 

PL000670 

PL000680 

PL000690 

PL000700 

PL0007I0 

PL000720 
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C 


CALL  CX41  (4107) 

CALL  PAGE  (8.5,11.) 

CALL  AREA2D(6.,6.) 

CALL  HEADIN  ('AMBIENT  NOISE  PLOT* , 20. 1 . , 1 ) 

CALL  POLAR( .01745.DR.3..3. ) 

CALL  POLORG  (MIN) 

CALL  GRIO(-30.2) 

CALL  MARKER  (15) 

CALL  CURVE(Y,Z. 360.1) 

CALL  REALNO(MIN,>0.2.8.2.7) 

CALL  REALNO(DRl,-fO,3.3,2.7) 

CALL  REALNO(DR2,+0.3.8,2.7) 

CALL  REALH0(DR3,^0,4.3.2.7) 

CALL  REALN0(DR4.«0.4.8.2.7) 

CALL  REALNO(DR5, ■'^0,5. 3.2.7) 

CALL  REALNa(MAX,-fO,5.8.2.7) 

CALL  MESSAG( 'NOISE  LEVELS  GIVEN  IN  DB  REF  lUPA' , 34. 1 . 0, -1 . 0 ) 
CALL  MESSAGCN',  1.3. 0.6.1) 

CALL  MESSAG('180'.3.2.8,-.2) 

CALL  MESSAG('090'.3.6.2.3.0) 

CALL  MESSAG( •270',3.-.5.3.0) 

CALL  ENDPL  (0) 

CALL  DOHEPL 

STOP 

END 


PL000730 

PL000740 

PL0007SO 

PL000760 

PL000770 

PL000780 

PL000790 

PL000800 

PL000810 

PL000820 

PL000830 

PLO00S40 

PL0008SO 

PL000860 

PL000870 

PL000880 

PL000890 

PL000900 

PL000910 

PL000920 

PL000930 

PL000940 

PL000950 

PL000960 

PL000970 
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FILE.  PL0T2 


FORTRAN  A1 


DAVID  C.  MCDONNELL,  LT  USN 

THIS  PROGRAM  USES  THE  DISSPLA  ROUTINE  TO  PLOT  THE  AG  DATA 
GENERATED  FROM  THE  OPTIMAL  PROGRAM.  THE  PROGRAM  READS  THE 
DATA  FROM  THE  FILE  AND  PLOTS  THE  DATA  ON  A  POLAR  PLOT. 

THE  SHIPS  HEADINGS  AREREPRESENTED  IN  DEGREES  TRUE  AND  THE 
AMPLITUDE  OF  THE  AG  IS  REPRESENTED  IN  OB  REF  1  UPA. 

VARIABLE  DECLARATION  SECTION 

INTEGER  l.J.K 

REAL  X(360),Y(360),Z(360),A.D,MIN,MAX,DR,DR1,DR2,ORS,TGT 
REAL  DR<i.DR5,STEP.Tl 
CHARACTERXZO  FILE$,FLN$ 

OBTAIN  THE  FILE  NAME  FROM  THE  USER  MHERE  THE  DATA  IS  STORED 

MRITECK,*) 'ENTER  FILENAME  AND  FILETYPE  HHERE  DATA  IS  STORED* 
READ(«, '(AZO)')  FILES 

A  MUST  BE  APPENDED  TO  THE  FRONT  OF  THE  FILENAME  FOR  THE 
DISSPLA  PROGRAM  TO  RECOGNIZE  IT. 

FLN$='/'//FILES//'A' 

INITIALIZE  VARIABLES 

A  =  0. 

MIN=100. 

MAX=0. 


PLOOOOIO 

PL000020 

PL000030 

PLOOOOAO 

PL000050 

PL000060 

PL0Q0070 

PL000080 

PL000090 

PLOOOlOO 

PLOOOllO 

PL000120 

PL000130 

PL000140 

PL000150 

PL000160 

PL000170 

PL000180 

PL000190 

PL000200 

PL000210 

PL000220 

PL000230 

PL000240 

PL000250 

PL000260 

PL000270 

PL000280 

PL00C290 


THE  FOLLOMING  DO  LOOPS  READ  THE  DATA  FROM  THE  FILE.  THIS  IS  DONE 
THICE  DUE  TO  THE  RECIPROCAL  HEADINGS  HILL  HAVE  THE  SAME  AG. 

DO  20  1*1,2 

OPEN  <UNIT*1,FILE*FLN$, STATUS* 'OLD') 

READ(l,)<)TGT 

IF  (TGT.LE.90)  Tl*90-TGT 
IF  (TGT.GT.90)  T1*450-TGT 
DO  10  J*l,180 
AsA+1 

READ  (1,)()  X(INT(A)} 

DETERMINE  THE  MINIMUM  AND  MAXIMUM  VALUE  OF  DATA  FOR  PLOT  PURPOSES 

IF  (X(INT(A)) .LT.MIN)  MIN*X( IHT( A) ) 

IF  (XdNTCA))  .GT.MAX)  MAX=X(INT(A) ) 

THE  DISSPLA  POLAR  COMMAND  PLOTS  ANGLES  IN  RADIANS  FROM  EAST. 
THEREFORE  DEGREES  TRUE  ARE  CONVERTED  TO  RADIANS  FROM  EAST  IN 
THE  FOLLOHING  IF  STATEMENTS. 

IF  (A.GT.0.AMD.A.LT.90)  THEM 
B=90-A 

ELSE  IF  (A.GE.90)  THEN 
B*450-A 
END  IF 

Y(INT(B))*REAL(B) 

Z<INT(B})*X(INT(A)) 

CONTINUE 
CLOSE  (1) 

CONTINUE 

SET  MIN/MAX  EQUAL  TO  +/-  10  FOR  PLOTTING  PURPOSES 

MIN  *  MIN  -  10. 

MAX  =  MAX  +  10. 

DR=INT((MAX-MIN)/6) 

STEP*2«DR 

IF  ((0-MIN).GE.(3*DR))  THEN 
DR4*0 

0R3*DR4-DR 


PLOOOSOO 

PLOOOSIO 

PLO00S20 

PL000330 

PLO00S40 

PL000350 

PL000360 

PL000570 

PL000380 

PL000390 

PL000400 

PL0004I0 

PL000420 

PL000430 

PL000440 

PL000450 

PL000460 

PL000470 

PL000480 

PL000490 

PLOOOSOO 

PLOOOSIO 

PL000S20 

PLO00S30 

PL000S40 

PLOOOS50 

PLO00S60 

PLO00S70 

PLO00S80 

PLO00S90 

PL000600 

PL000610 

PL000620 

PL000630 

PL000640 

PLO006S0 

PL000660 

PL000670 

PL000680 

PL000690 

PL000700 

PL000710 

PL000720 
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0R2-0R3-DR 
DR1=DR2-DR 
MIN-DRl-DR 
DR5=DR4^DR 
MAX=DR5^DR 

ELSE  IF  ((0-MIN).LT.(3XDR))  THEN 
DR2  =  0 

DR1=DR2-0R 
MIN=DR1-DR 
DR3=DR2+DR 
0R4=0R3«DR 
DR5=DR4+DR 
MAX=DR5>0R 
END  IF 

THE  PLOTTING  ROUTINE  USES  DISSPLA  COMMANDS  TO  PLOT  THE  DATA 
IN  DEGREES  TRUE  HITH  LEVELS  GIVEN  IN  DB  REF  lUPA. 

CALL  SHERPAC  •MCT5DATAS*A',3) 

C  CALL  CX41(4107) 

CALL  PAGE  (8.5,11.) 

CALL  AREA2D(6..6.) 

CALL  HEADIN  ('SNR  GAIN  RELATIVE  TO  BROADSIDE  HEADING* , 38, 1 I ) 
CALL  POLAR(.0174S,STEP,3.,3.) 

CALL  POLORG  (MIN) 

CALL  GRID  (-30,2) 

CALL  MARKER(4) 

CALL  CURVE(T1. MAX, 1,-1) 

CALL  MARKER  (15) 

CALL  CURVE  (Y, 2,360,1) 

CALL  REALNO(MIN,-»0,2.8,2.7) 

CALL  REALNO(OR1,-I-0,3.3,2.7) 

CALL  REALN0(DR2,^0,3.8,2.7) 

CALL  REALN0(DR3,>0,'4.3,2.7) 

CALL  REALN0(DR4,-f0,4.8,2.7) 

CALL  REALNO(DR5,-»0,5.3,2.7) 

CALL  REALMQ(MAX,+0,5.8,2.7) 

CALL  MESSAG( 'LEVELS  GIVEN  IN  DB  REF  lUPA* ,28, I .3,-1 . 0) 

CALL  MESSAG( 'N', 1,3. 0,6.1) 

CALL  ME5SAG('180',3,2.8,-.2) 

CALL  MESSAG('090*,3,6.2,3.0) 

CALL  MESSAG('270',3,-.5,3.0) 

CALL  MESSAGCX  =  TARGET  TRUE  BEARING* , 23, 1 . 3, -1 . 4) 

CALL  REALN0(TGT,+0,4.2,-1.4) 

CALL  ENDPL  (0) 

CALL  DONEPL 

STOP 

END 


PL000730 

PL000740 

PL000750 

PL000760 

PL000770 

PL000780 

PL000790 

PL000800 

PL000810 

PL000820 

PL000830 

PL000840 

PLO0Q850 

PL000860 

PL000870 

PL000880 

PL000890 

PL000900 

PL000910 

PL000920 

PL000930 

PL000940 

PL000950 

PL000960 

PL000970 

PL000980 

PL000990 

PLOOIOOO 

PLOOlOlO 

PL001020 

PL001030 

PL001040 

PL001050 

PL001060 

PL001070 

PL001080 

PL001090 

PLOOllOO 

PLOOlllO 

PL001120 

PL001130 

PL001140 

PL001150 

PL001160 

PL001170 

PL001180 

PL001190 

PL001200 
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APPENDIX  B 
ADDITIONAL  EXAMPLES 


Examples  of  the  "Optimal"  program,  along  with  a  discussion  of 
each,  are  included  in  this  appendix.  The  first  three  examples  are  test 
cases  which  partially  validate  the  program,  and  the  other  examples 
show  the  applicability  of  the  program,  depicting  the  SNR  gain 
achievable  in  some  cases  when  optimal  heading  selection  is  used. 

A.  EXAMPLE  1 

Example  1  uses  an  omnidirectional  ambient  noise  field  which 
has  a  level  of  50  dB  as  seen  in  Figure  9.  In  this  case  the  broadside 
heading  should  be  the  optimal  heading  due  to  the  narrow  beamwidth 
of  the  broadside  beams.  Figure  10  shows  the  results  of  the  "Optimal" 
program  for  this  case  with  the  suspected  target  location  bearing  090. 
Notice  the  highest  SNR  gain  is  for  the  headings  oriented  towards  the 
north  and  the  south,  the  broadside  headings.  As  the  linear  array's 
orientation  is  changed  away  from  the  broadside  headings,  the  SNR 
gain  becomes  negative.  Thus,  the  array's  performance  is  degraded  by 
up  to  about  4  dB  when  the  array  is  oriented  to  place  the  target  in  the 
endfire  beams.  When  the  array  is  oriented  to  place  the  target  just 
off  of  endfire,  the  array's  performance  is  even  worse  than  that  of  the 
endfire  headings.  This  is  due  to  the  fact  that  the  main  lobe  and  the 
ambiguous  beam  overlap  and  a  large  heart  shaped  beam  pattern  is 
formed  which  has  a  larger  beamwidth  than  that  of  the  endfire  beam. 
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AMBIENT  NOISE  PLOT 


N 


180 

NOISE  LEVELS  GIVEN  IN  DB  REF  l^JPfl 


Figure  9.  Omnidirectional  ambient  noise  field  at  150  Hz.  Levels  are 
given  in  dB  ref  1  jiPa. 
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SNR  GAIN  RELATIVE  TO  BROADSIDE  HERDING 


N 


180 

LEVELS  GIVEN  IN  DB  REF  IJJPR 
X  -  TARGET  TRUE  BEARING  90. 


Figure  10.  SNR  gain  for  all  true  headings  relative  to  the  broadside 
heading. 
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B.  Example  2 

Example  2  provides  another  test  case  which  is  presented  to 
partially  validate  the  "Optimal”  program.  An  ambient  noise  field 
which  has  a  very  high  noise  level  at  just  one  true  bearing,  090,  was 
created  and  is  displayed  in  Figure  11.  With  the  suspected  target 
location  to  the  west,  the  broadside  headings  of  000  and  180  would 
place  the  ambiguous  beam  in  the  direction  of  the  high  noise  level  at 
090  and  the  array  performance  would  be  severely  degraded.  Any 
heading  which  does  not  place  the  array's  ambiguous  beam  in  the 
direction  of  the  high  noise  level  will  have  a  significant  SNR  gain 
relative  to  the  broadside  heading.  As  can  be  seen  in  Figure  12,  the 
SNR  gain  for  almost  all  headings  is  anywhere  from  11  to  18  dB 
higher  than  the  SNR  for  the  broadside  heading.  Figure  12  shows  the 
results  of  the  array's  performance  as  the  ambiguous  beam  passes 
through  the  high  noise  level  located  at  090.  Note  the  curved  pattern 
of  the  SNR  gain  plot  in  Figure  12  showing  the  inverted  beam  pattern 
of  the  array  for  the  headings  which  place  the  ambiguous  beam  in  the 
high  noise  level  area.  Additionally,  total  noise  level  calculations 
using  Eq.  (24)  verify  the  program's  result  of  the  18  dB  SNR  gain  for 
the  headings  just  off  broadside. 


AMBIENT  NOISE  PLOT 


NOISE  LEVELS  GIVEN  IN  OB  REF  l^Pfl 


Figure  11.  Highly  directional  ambient  noise  field  at  150  Hz.  Levels 
are  given  in  dB  ref  1  jiPa.  Note  the  high  noise  level  located  at  a  true 
bearing  of  090. 


40 


SNR  GRIN  RELATIVE  TO  BROADSIDE  HERDING 


LEVELS  GIVEN  IN  DB  REF  i;jPR 
X  -  TARGET  TRUE  BEARING  270. 


Figure  12.  SNR  gain  for  all  true  headings  relative  to  the  broadside 
heading  in  a  highly  directional  noise  field.  Note  the  inverted  beam 
pattern  for  array  orientation  of  000  and  180. 
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C  EXAMPLE  3 

Example  3  is  another  highly  directional  ambient  noise  field 
which  is  used  to  partially  validate  the  "Optimal"  program.  In  this 
case  very  high  noise  levels  of  80  dB,  compared  with  the  rest  of  the 
noise  field,  are  located  between  true  bearing  of  160  and  200  as  seen 
in  Figure  13.  With  the  suspected  target  location  towards  the  north, 
the  broadside  headings  of  the  array  place  the  ambiguous  beam  in  the 
direction  of  the  high  noise  levels.  The  array's  performance  will  be 
degraded  for  headings  which  place  the  ambiguous  beam  in  the  high 
noise  level  area.  Figure  14  shows  the  results  of  the  "Optimal" 
program  for  this  case.  The  SNR  gain  for  all  headings  which  do  not 
place  the  ambiguous  beam  in  the  high  noise  level  area  perform  much 
better  than  the  broadside  headings  of  090  and  270.  Again  the 
inverted  beam  pattern  of  the  array  is  represented  centered  about 
the  broadside  headings.  The  constant  SNR  level  of  0  dB  gain  results 
when  the  entire  ambiguous  beam  is  in  the  high  noise  level  area.  As 
the  ambiguous  beam  starts  to  exit  out  of  the  high  noise  level  area, 
the  SNR  gain  gradually  increases  up  to  27  dB,  at  which  point  none  of 
the  ambiguous  beam  is  in  the  high  noise  level  area.  Total  noise  level 
calculations  using  Eq.  (24)  verify  the  program's  result  of  a  27  dB  gain 
for  the  headings  just  off  broadside,  which  do  not  place  the  ambiguous 
beam  in  the  high  noise  level  area. 
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AMBIENT  NOISE  PLOT 


Figure  13.  Highly  directional  ambient  noise  field  at  150  Hz.  Levels 
are  given  in  dB  ref  1  ^Pa.  Note  the  high  noise  level  located  towards 
the  south. 
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SNR  GRIN  RELATIVE  TO  BROADSIDE  HERDING 


LEVELS  GIVEN  IN  DB  REE  lJUPR 
X  -  TARGET  TRUE  BEARING  360. 


Figure  14.  SNR  gain  for  all  true  headings  relative  to  the  broadside 
heading  in  a  highly  directional  noise  field. 
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D.  EXAMPLE  4 

Figure  15  displays  an  actual  ambient  noise  field.  Noise  levels 
vary  from  a  high  of  about  60  dB  towards  the  east-southeast  to  a  low 
of  about  45  dB  towards  the  west.  The  "Optimal"  program  was  used  to 
analyze  for  a  suspected  target  location  towards  bearing  315.  As  seen 
in  Figure  16,  gains  of  up  to  3  dB  can  be  obtained  by  choosing  the 
appropriate  array  headings  in  this  case.  For  instance,  choosing  a 
array  heading  of  about  285  will  produce  a  SNR  gain  of  3  dB  relative 
to  the  broadside  heading  of  045. 
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AMBIENT  NOISE  PLOT 


N 


180 

NOISE  LEVELS  GIVEN  IN  OB  REF  IjJPfl 


Figure  15.  An  actual  directional  ambient  noise  field  at  150  Hz. 
Levels  are  given  in  dB  ref  1  fxPa. 
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SNR  GRIN  RELATIVE  TO  BROADSIDE  HERDING 

N 


180 


LEVELS  GIVEN  IN  DB  REF  IJJPR 
X  -  TARGET  TRUE  BEARING  315. 


Figure  16.  SNR  gain  for  all  true  headings  relative  to  the  broadside 
heading  in  an  actual  directional  noise  field. 


E.  EXAMPLE  5 

Figure  17  is  an  example  of  an  actual  ambient  noise  field  which 
has  two  high  noise  level  peaks  of  60  dB  at  true  bearings  of  about  050 
and  085.  Additionally,  the  noise  levels  toward  the  west-southwest 
are  lower  with  levels  of  only  about  45  dB.  The  "Optimal"  program 
was  run  for  this  case  and  the  results  are  displayed  in  Figure  18.  Note 
the  large  variation  in  the  SNR  gain  with  relatively  small  changes  in 
array  heading.  For  instance,  for  an  array  heading  of  330,  a  4  dB  SNR 
gain  is  achieved,  but  for  an  array  heading  of  about  335,  the  "'NR  gain 
is  -2  dB.  A  6  dB  difference  in  array  performance  is  obtained  with 
only  a  change  of  5  degrees  in  array  orientation.  The  two  peaks  in  the 
ambient  noise  field  are  responsible  for  the  high  variability  in  the  SNR  ’ 
gain  in  this  case. 
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AMBIENT  NOISE  PLOT 


N 


180 

NOISE  LEVELS  GIVEN  IN  DB  REF  IJJPR 


Figure  17.  An  actual  directional  ambient  noise  field  at  150  Hz. 
Levels  are  given  in  dB  ref  1  uPa. 
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SNR  GAIN  RELATIVE  TO  BROADSIDE  HEADING 


LEVELS  GIVEN  IN  DB  REF  IJUPA 
X  -  TARGET  TRUE  BEARING  270. 


Figure  18.  SNR  gain  for  all  true  headings  relative  to  the  broadside 
heading  in  an  actual  directioncil  noise  field. 
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